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Abstract 

We study the fluctuations of the interface, in the steady state, of the Surface Relaxation Model 
(SRM) in two scale-free interacting networks where a fraction q of nodes in both networks interact 
one to one through external connections. We find that as q increases the fluctuations on both 
networks decrease and thus the synchronization reaches an improvement of nearly 40% when q = 1. 
The decrease of the fluctuations on both networks is due mainly to the diffusion through external 
connections which allows to reducing the load in nodes by sending their excess mostly to low- 
degree nodes, which we report have the lowest heights. This effect enhances the matching of 
the heights of low-and high-degree nodes as q increases reducing the fluctuations. This effect is 
almost independent of the degree distribution of the networks which means that the interconnection 
governs the behavior of the process over its topology. 
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I. INTRODUCTION 


In the last decades the study of complex networks has been growing strongly due to 
the large number of systems that exhibit this type of structures. A complex network is a 
set of nodes that are connected by internal links and the most fundamental property that 
characterizes its topology is the degree distribution P(fc), which represents the probability 
that a node has k neighbors or connectivity k. It has been found that many real systems 
such as social, communication and biological networks present a degree distribution given 
by P{k) ~ k~^, where A is the exponent of the power law and kmin <k< k^ax, where kmin 
and kmax are the minimum and maximum degree of the network. These kind of networks 
are called Scale Free (SF) and one of its most important features is that are in general very 
heterogeneous, i.e. most nodes of the network have a low connectivity while only a few 
have a high connectivity (hubs). In recent years the study of synchronization processes in 
isolate complex networks has been increasing because of its importance in neurobiology 
and population dynamics p, 7|. A common theoretical approach to study synchronization in 
complex networks is to map this process onto an interface growth model by assigning to each 
node a scalar held hi, with i = 1, ■ ■ ■, N, where N is the size of the network. This scalar held 
could represent, for example, the amount of load on a node in the problem of distributed 
parallel computing on processors. Without loss of generality we will relate the scalar held 
to a set of heights on the interface. The most relevant magnitude that characterizes the 
interface is W{t) = W, which represents the huctuations of the scalar held around its mean 
value on the network, given by 


N 


1/2 
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where {h) = A hi is the average of the scalar held over the nodes at time t and { } 
is the average over diherent network realizations. The roughness evolves in time until it 
saturates at the steady state. In the saturation regime Wg is a constant that depends only 
on A. In complex networks the synchronization of the system is related to the roughness in 
the saturation regime js 


16|. 


One of the most simple and used mode 
is the Surface Relaxation Model (SRM) [6 
model largely used in euclidean interfaces 


s to study synchronization in complex networks 


12] which is a modication of the classic Family 
ITj j. In this model, at each time step, a node i 
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is randomly chosen and the node with the lowest height between the chosen node and all 
its neighbors evolves increa sing its height. It has been found that for isolated SF networks, 


with A < 3, Ws ~ IniV 8, 10, 18| until a critical value N = N*, af 


independent of iV, which means that the system becomes scalable 


er which W. becomes 


18l |. Although it was 


an interesting result, many real systems are not isolated but interacting with other systems 
instead. This means that a process that develops in one network can be affected by a process 


developing in another and vice versa 
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35| . Such is the case of epidemic models where 


the interaction between networks make it very harmful for the healthy populations because 
the interaction increases the theoretical risk of infection compared with the same process in 
isolated networks 


36 




These interacting systems can be modeled as networks that interact through external 
links that connect nodes that belongs to different networks. Now the question is, does 
synchronization in interacting networks performs worse or better than synchronization in 
isolated networks?. In order to answer this question, in this letter we study the synchroniza¬ 
tion of two SF networks with the same size N that interact through a fraction q of nodes 
connected, one by one, between them. The model used is the SRM model which we adapt to 
interacting networks and study the effect of the interaction parameter q on the fluctuations 
in both networks. 


II. MODEL 

In our model, two uncorrelated SF networks, called A and B, with the same size N and 
exponent and \b respectively are built using the Molloy-Reed algorithm ^ disallowing 
self loops and multiple connections. As we need a single interface on each network, in order 
to ensure that we have a single component we use kmin = 2 [43|. Each node z G a, with 
a = A,B, has a connectivity kf and we denote the set of its neighbors by vf. In order 
to build the external connections between the networks we connect by simplicity, the hrst 
q N nodes in A one by one with the hrst q N in R, where q is the interaction parameter 
with 0 < g < 1. Notice that if both networks are uncorrelated this procedure is the same 
as connecting a fraction q of nodes at random. We dehne the vector M, where Mj, with 
i = 1, is equal to 1 if the node i in A has an external connection with i in R, and 

Mj = 0 otherwise. In order to simplify the growth rules we choose random initial conditions 
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FIG. 1: Schematic of the rules of the model in a one dimensional Euclidean lattice. The particles P 
are only dropped on the network A in this scheme. The numbers represent different nodes in each 
network and the red dotted lines represent the external connections between nodes of different 
networks. In this case = 11 and q = 3/11. The arrows indicate the path that the particles 
follow, which goes from the node where the particle was originally dropped, to the node where 
the particle gets finally deposited. The height of the nodes is measured from the upper line of the 
boxes that represent the numbers assigned to the nodes in each network. 

for the scalar field in the interval [0,1], hence, we avoid the cases in which different nodes 
have eqnal heights, as we are only interested in the saturation regime on both networks 
where the initial condition plays no role. 

The evolution rules of the interface growth are given as follows: 

1. A network a (with a = A, B) is chosen with probability 1/2 and then a “particle”, 
which represents the load, is dropped in a node i selected randomly in a. 

2. The particle diffuses to the node e that is the node with the lowest height between the 
node i and its neighbors n“. 
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FIG. 2: as a function of q for A (Q) and B (□) with = 2.6 and = 3. In the inset 

Wg/Wg{q = 0) as a function of q. 

3. If Mg = 0 or Mg = 1 and ^ a < & (3 (with (3 ^ a) the particle is deposited in 

e G a. Otherwise the particle diffuses to the network (3 and is deposited in the node 
with the lowest height between e and its neighbors nf. 

Thus if we denote ^ & a a.s the node where the particle is hnally deposited, then h" = 
h" + 1. At each Monte-Carlo step the time is increased by \/2N. In Fig. [U we show a 
schematic of the rules of the process for the case of a one-dimensional lattice. 


III. RESULTS AND DISCUSSIONS 


We are interested in the behavior of the fluctuations in the steady state of both networks 
with the system size above N* 1^, value for which the system is scalable. For isolated SF 
networks this regime for A < 3 is close to N* 10^ [^. We check that the nature of this 

regime is due to the distribution of internal connectivities and that it is almost not affected 
by the interaction parameter q. Thus in our research we use N = Na = Nb = 3 x 10^ 
in order to ensure that we are in the scalable regime. We will show our results only for 
Aa = 2.6 and A^ = 3, because all the other combinations of the exponents A in 2.5 < A < 3 
give qualitatively the same results. We compute the fluctuations in the saturation regime 
of both networks IT" = \J jj ~ 'with a = A, B and in Fig. [2] we show Wg 

as a function of q. It is clear that as the interaction parameter q increases, Wg decreases 
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FIG. 3: — (/i“) as a function of A; for g = 0 in ^ (A) and B (y), q = 0.5 in A (x) and B (□) 

and q = 1 iw A (Q) and B (*). The inset is an amplification of /i^ — (/i“) for the lowest values of 
A: for g = 0 in A (A) and g = 1 in ^4 (Q) 

in both networks, which implies that the synchronization improves. From the plot we can 
also observe that as q increases, the difference between Wg and W§ becomes smaller, which 
means that the synchronization in each network becomes mainly controlled by q and not 
by the internal degree distributions. In the inset of Fig. |2] we show Wg /Wg{q = 0) as a 
function of q. It can be seen how is the rate of improvement in the synchronization with the 
increment of q. For example, for q = 0.3 the synchronization enhances approximately 20% 
and for g = 1 around 40%. It is worth pointing out that this rate decreases as q increases, 
and this means that the effect of the optimization gets less signihcant as the networks have 
more interconnections between them. 

To understand the effect of the interaction parameter q on the optimization of the process 
we compute the difference between the average height of nodes with degree k, denoted by 

and the mean value of the height (h“) as a function of k. In Fig. |3]we show — (h“) 
as a function of k for different values of q. We can see that, for g = 0 the heights of low 
connectivity nodes, which are the majority in SF networks, are closer to the average height 
of the network than the heights of high degree nodes, which are above (h"). This means 
that hubs are usually overloaded because all their neighbors send them their excess of load, 
affecting negatively the synchronization of the system. However, as the factor g increases 
the height of the hubs decreases, approaching to (h") and becoming independent of k for 
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FIG. 4: Rates at which a particle spreads from the network a to the network (3 as a. function of 
q. With uaa (O)) ^ab I'BB (x) and vba (A). From the definition it is clear that vaa+^ab = 1 
and vba + ^bb = 1 - 



FIG. 5: Degree distributions Pap{k) of the particles which get deposited on a node with connectivity 
k on the network /? after being dropped on network a for <7 = 0.5. With Paa (O)) Pab Pba 
(x) and Pbb (A). 

q = 1. In the inset of Fig. |3l we can see an amplification of the behavior of for the nodes 
with the lowest connectivity. We can see that as q increases their heights also approach to 
the average value. These results imply that as we increase the factor q hubs are no longer 
overloaded and therefore their excess of load is now absorbed by low connectivity nodes. 

In order to understand the diffusion process between networks we want to know how 
frequently the load spreads from one network to another and how it gets distributed after 
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TABLE I: Dispersion cJq,^ of the distributions for different values of q 


q 


= 2.6 As = 

= 3.0 

(^AA 


(^AB 

BA 

0 

22.1 

7.5 

- 

- 

0.1 

22.4 

7.7 

2.2 

3.3 

0.5 

23.9 

8.2 

2.7 

4.5 

1.0 

25.7 

8.6 

3.0 

5.7 


diffusion. Hence we measure the rates at which a particle dropped in the network a 
gets hnally deposited in the network (3. In Fig. |4] we plot Pap as a function of q. When 
a = (3 the rate always decreases with q, because the particles have more chances to cross 
to the other network. We can see that the rate is always much bigger when a = f3 than 
in the case of a 7^ (3, which means that the particles tend to diffuse in the same network 
most of the time. The fact that a small portion of load that crosses between networks is 
enough to enhance the synchronization is due to the fact that the dynamics of the SRM 
in both networks are coupled due to the interaction. Another observation is that all the 
rates do almost not depend on the degree distribution of each network. The majority of 
nodes with Mj 7^ 0 have an external connection with low connectivity nodes, due to the fact 
that these are selected at random and are the most common ones in SF networks. The fact 
that the amount of nodes with low connectivity does not change much with variations in 
Xa for 2.5 < Aq, < 3, and also because the diffusion between networks depends directly on 
the externally connected nodes whether a particle crosses to another network or remains in 
the same, make that the rates have almost no dependence on the exponents of the original 
degree distributions. 

How the load gets distributed after the diffusion process? In order to answer this question 
we compute the probability Pap{k), dehned as the probability that a particle dropped in 
the network a gets deposited in a node with degree k in the network (3. In Fig. [5] we 
plot Pap{k) with a,/? = A, H for g = 0.5. Also in Table [T] we report the dispersion aap of 
these distributions, which quantihes the heterogeneity of the deposition process, for different 
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FIG. 6: Rates and at which a particle spreads from a to network and gets directly 
deposited or deposited on a neighboring node respectively. With (0)j ^ab ^BA (^) 

^BA (A). 

values of g and for the isolated networks. From the plot we can see that the probabilities 
Paa{k) for a = A, B are very heterogeneous and have a similar dispersion to the dispersion 
in the isolated networks. This means that when a particle stays on the network where it was 
originally dropped hnds an environment with a similar heterogeneity than in the isolate SF 
network. Moreover the distributions Pag{k) for a ^ f3 are more homogeneous than in the 
case a = (3 and thus have a lower dispersion. This means that when a particle crosses from 
one network to another it hnds a more homogeneous neighborhood than in the case in which 
stays in the same network, and as we will show below low-degree nodes are hlled more often 
than high-degree nodes. Also we can see that the dispersion, for all cases, slightly grows 
with q, and this is due to the fact that hubs are less overloaded and can participate more 
often in the diffusion process. 

We want to understand the reason that makes the load get rather deposited on low- 
connectivity nodes when it crosses to a different network. In order to explain this effect we 
study the diffusion process when the load crosses to the other network. When a particle 
spreads from one network to another, it can be directly deposited on a node connected by 
the external connection or it can be deposited in one of its neighbors. To understand which 
of these two scenarios is more probable, we compute and which are the rates at 
which a particle that spreads from one network to another gets deposited directly {d) or in 
a neighboring (n) node respectively. Notice that = 1. In Fig. [6] we plot these 
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FIG. 7: Probabilities P^fsik) and that after an external diffusion, the particle gets deposited 

directly or gets deposited on a neighboring node with degree k respectively for q = 0.5 . With Pj^^ 
{0),P1b (□), F^^(x)andP^^ (A). 

rates as a function of q and we can see that is more important than which means 
that for any g, most of the times the particles that cross from one network to another get 
directly deposited. 

In order to explain the last observation, we dehne P^isik) and P^qik) as the probabilities 
that, after an external diffusion, the particle gets directly deposited or gets deposited on 
a neighboring node with degree k respectively. In Fig. [7] we plot these probabilities for 
q = 0.5. We can see that the probabilities P^isik) are more homogeneous than P^pik) and 
that nodes with low degree are the ones which receive the majority of the particles that are 
directly deposited. We can also see that the probabilities P^isik), which contemplate the 
scenario of neighboring deposition, have a wider spectrum, which agrees with the fact that 
in SF networks there are a few nodes with high degrees that receive the load by diffusion 
from their neighbors. This mechanism reduces the fluctuations, due to a matching of the 
heights of low-degree and high-degree nodes that is more efficient as q increases. Finally the 
system is optimally synchronized for q = 1. 


IV. CONCLUSIONS 

We study the synchronization in two SF networks where the dynamic of growth is ruled 
by a modihed SRM model. We study the fluctuations in the steady state of each network 
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as a function of the interacting parameter q and we find that the synchronization of each 
network improves when the interconnection between them increases. This improvement in 
both networks is about a 40% better than in isolated networks when q = 1, which is an 
important value regarding the decrease of the fluctuations. However, we show that the rate 
of this improvement decrease with the number of interconnections and this can be an useful 
result in future research to determine if a larger interconnection between networks apart 
from g = 1, is worth it. 

On the other hand, the improvement in the synchronization is due mainly to the diffusion 
through external connections. We also found that the majority of the particles that travel 
through the external connections are directly deposited in nodes with low connectivity, which 
are the majority in SF networks. We observe that in average these nodes usually have a the 
lowest heights in a SF network. Then when q increases, the height of the nodes with low 
connectivity increases compared to the mean value and the difference with the height of the 
hubs decreases. We also found various distinctive characteristics of the model, such as the 
fact that the synchronization and the heights of the hubs become almost independent of the 
internal degree distribution as q increase, or the fact that the percentage of particles that 
diffuse between networks only depends on q and not on the internal degree distribution of 
the networks. This result indicates that for high values of q the behavior of the system is 
ruled by the interconnection between networks and not by the topology of the system. 
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